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The present paper proposes a novel interpretation of the widely scattered states (called synchro- 
nized traffic) stimulated by Kerner's hypotheses about the existence of a multitude of metastable 
states in the fundamental diagram. Using single vehicle data collected at the German highway 
Al, temporal velocity patterns have been analyzed to show a collection of certain fragments with 
approximately constant velocities and sharp jumps between them. The particular velocity values in 
these fragments vary in a wide range. In contrast, the flow rate is more or less constant because its 
fluctuations are mainly due to the discreteness of traffic flow. 

Subsequently, we develop a model for synchronized traffic that can explain these characteristics. 
Following previous work (I. A. Lubashevsky, R. Mahnke, Phys. Rev. E 62, 6082, 2000) the vehicle 
flow is specified by car density, mean velocity, and additional order parameters h and a that are 
due to the many-particle effects of the vehicle interaction. The parameter h describes the multilane 
correlations in the vehicle motion. Together with the car density it determines directly the mean 
velocity. The parameter a, in contrast, controls the evolution of h only. The model assumes that 
a fluctuates randomly around the value corresponding to the car configuration optimal for lane 
changing. When it deviates from this value the lane change is depressed for all cars forming a 
local cluster. Since exactly the overtaking manoeuvres of these cars cause the order parameter a to 
vary, the evolution of the car arrangement becomes frozen for a certain time. In other words, the 
evolution equations form certain dynamical traps responsible for the long-time correlations in the 
synchronized mode. 



I. ANOMALOUS PROPERTY OF THE 
SYNCHRONIZED MODE 



Although the motion of individual vehicles is controlled 
by the motivated driver behavior rather than by sim- 
ple physical laws the car ensembles on highways exhibit 
phenomena widely met in physical systems. Namely, 
the existence of various states (self-sustained steady-state 
modes of traffic flow), their coexistence, the phase tran- 
sitions, etc. (reviewed in [[jj g]), actually forms a novel 
branch of physics. Following the pioneering works by 
Lighthill & Whitham Q], and Richards || the state of 
traffic flow on a highway is specified by the vehicle den- 
sity p and the mean velocity v or, equivalently, by the 
traffic flow rate q = pv. In other words, the vehicle den- 
sity p and the flow rate q are regarded as the complete set 
of state variables, leading to a correspondence between 
a point in the pg-plane with a particular state of traf- 
fic flow. The main assumption adopted in the classical 
theories (reviews can be found in |(| fj], ||, |[ |l(J) is the 
existence of a relationship between the mean velocity and 
vehicle density, v = $(p). So, the curve q = p${p) in the 



pq-planc is called fundamental diagram which character- 
izes the possible states of traffic flow. 

According to traffic flow data analyzed recently by 
Kerner & Rehborn jyj |§ Kerner @, |]|, and also 
Neubert et al. |l(| there are three distinctive states of 
multilane traffic flow, the free flow (F), the synchronized 
mode (S), and wide moving jams (J). This endows the 
multilane traffic flow with a variety of properties. In par- 
ticular, it has been demonstrated that the self- formation 
of moving jams mainly proceeds via the sequence of two 
phase transitions F — ► S — ► J p4| . Both of them are 
of the first order, i.e. exhibiting breakdown, hysteresis, 
and nucleation effects |ll| . The present paper focuses its 
attention on a unique property of the synchronized mode 
itself, its complexity |l2| . In contrast to the free flow, 
the synchronized mode matches a two-dimensional do- 
main on the pg-plane rather than a curve on it (Fig. ^). 
It means that a track made up of the empirical points 
{p(t),q(t)} obtained at sequential time moments unpre- 
dictably fills a two-dimensional domain. In this regard 
the synchronized mode is also referred to as a widely 
scattered traffic state for which a fundamental diagram 
in the form of a one-dimensional curve did not exist. 
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FIG. 1: Illustration of the experimental data by ]12| showing 
possible states of the free flow and the synchronized mode in 
the pq-plane. 



The synchronized mode is characterized by strong cor- 
relations in the vehicle motion at different lanes ordering 
the car arrangement across the road. To describe this 
state a model has to take into account particular details 
of the multilane car interaction. Whence it follows that 
the vehicle density p and the mean velocity v averaged 
over all the lanes across the highway form an incomplete 
set of state variables. One of the possible ways to tackle 
this problem is to ascribe to each lane its own vehicle den- 
sity and mean velocity. Several models for multilane traf- 
fic flow applying to such an approach have been proposed, 
in particular, gas-kinetic theo ry Jlj] , jig) , [l9| , |2Q , g]J , |||] 
and compressible fluid model |23| . For a more detailed 
survey, including also the cellular automata multilane 
models, we refer to || [?J ||, [| [l7J]. Here we only note 
that in all such models the cause of the free flow insta- 
bility is related to a delay in the driver's response to 
variations of the headway distance. This instability gives 
rise to the traffic flow pinch leading to the formation of 
a new congested phase of vehicle motion. 

In spite of the variety of available models for multi- 
lane traffic flow up to now no generally accepted expla- 
nation for the widely scattered states of the synchronized 
mode has been proposed. There are various points of 
view on the mechanism responsible for this phenomenon. 
In particular, it is a mixture of different vehicle types 
(cars and trucks) a heterogeneity in the head- 

way distance [^6) as well as in the highway structure |27| , 
changes in the behavior of "frustrated" drivers |2£| , anti- 
cipation effects [^9[ |o[ 31, |3^| , non-unique equilibrium 
solutions of the Prigoginc-Hcrman kinetic equation [p3[, 



the existence of certain plateaus in the dependence of the 
optimal velocity on the headway distance |54], and 
the variety of oscillating metastable states with different 
wavelengths |36, 37 . The requirement of multilane mod- 
els to describe synchronized traffic flow has already been 
stated in |lj| [38| using hydrodynamic as well as kinetic 
equations. 

Commenting on these models we would like to say the 
following. First, most of them apply to a one-dimensional 



representation of traffic flow, discrete or continuum one. 
Therefore they do not take into account the strong mul- 
tilane vehicle interaction which is an essential feature of 
the synchronized mode. Such a theory can explain the 
complex behavior of traffic flow on a single-lane road 
rather than on multilane highways or can describe the 
congested multilane traffic flow providing that the multi- 
lane synchronization is a minor effect. Second, according 
to the analysis of single- vehicle data in jl(| the synchro- 
nized mode is mainly singled out by strong multilane 
correlations in the vehicle velocity rather than the for- 
mal equality of speed at which cars moving on different 
lanes. In particular, in the observed free flow regime the 
averaged velocities at various lanes differ from one an- 
other only slightly. Therefore the approach based on the 
description of multilane traffic flow in terms of individ- 
ual streams along different lanes with their own mean 
velocities and vehicle densities seems to be doubtful. 

Following the classification proposed by Kerner & 
Rehborn |12] there are three distinctive kinds of syn- 
chronized flow: (i) stationary and homogeneous states 
where both the average speed and flow rate are ap- 
proximately constant during a fairly long time interval, 
(ii) states where only the average vehicle speed is sta- 
tionary ('homogeneous-in-speed-states'), and (Hi) non- 
stationary and non-homogeneous states. It should be 
noted that on macroscopic scales congested traffic ex- 
hibits a large variety of different phenomena, causing a 
more detailed classification |3^, [ll], [l2| (see also re- 
view ||). Recently, Kerner formulated several hypothe- 
ses about the synchronized mode |], |[ ||, [l2| to ex- 
plain the properties of the widely scattered states. In 
particular, regarding the synchronized traffic of types 
(i) and (ii) he assumed it to contain a whole multitude 
(continuum) of possible congested traffic states |l^, [l3| , 
which are stable with respect to infinitesimal perturba- 
tions. This leads to continuous spatial-temporal transi- 
tions between these states p], p[ (43|. However, sim- 
ulation models and further empirical evidence for these 
hypotheses are still to be found |J . 

According to the prevailing notion (see, e.g. Ref. |J) the 
synchronized mode is mainly identified by the absence of 
a direct relationship between the locally averaged vehicle 
velocity and density provided the case of the "stop-and- 
go wave" pattern has been excluded from consideration. 
Therefore, the main part of the theoretical models re- 
lates the synchronized mode to the instability of the ho- 
mogeneous traffic flow caused by the driver delay in the 
response to variations in the headway distance. These 
models explain the existence of the widely scattered traf- 
fic states taking into account the complex behavior of the 
developed spatial structure. From our point of view this 
is a quite natural approach that describes only the syn- 
chronized mode of type (Hi), called "oscillating congested 
traffic" ||. However, to justify Kerner's hypotheses con- 
cerning types (i) and (ii) in the framework of such an 
approach a traffic flow model has to admit a large vari- 
ety of stable nonhomogeneous solutions that either are 
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FIG. 2: The analyzed section of German highway Al and position of the loop detectors Dl - D4. Arrows indicate the direction 
of vehicle motion. 



fixed in the space or move as a whole with a constant 
speed. Such a model has been proposed in paper |3(| 
based on the car-following approach where the existence 
of the variety of stable nonhomogeneous solutions dif- 
fering in wavelength has been found numerically. It is 
an open question whether this variety is due to the pe- 
riodic boundary conditions adopted in the model or is 
an inherent property. At least under random noise, the 
model demonstrates a certain selection of a single stable 
solution, so, it seems to describe actually type (Hi) of 
the synchronized mode or a certain mixture of types (ii) 
and (in). We would like to underline once more that the 
cause of such a congested traffic state is related to the de- 
lay in the driver response giving rise to the instability of 
homogeneous traffic flow and the multilanc synchroniza- 
tion is a secondary effect only. Besides we note paper ]37j 
that models the widely scattered state region by apply- 
ing actually the idea used also by Kerner p|, p|, 3j, [43| in 
justifying his hypotheses. Namely it is assumed that in 
congested traffic drivers try to maintain a constant time 
to collision. The latter is possible within a wide interval 
of headway distances if the velocity difference between 
cars is sufficiently small. However this effect is a conse- 
quence of a certain degeneration of the adopted model. 
Introducing other terms they can again produce the os- 
cillating behavior of the developed nonhomogeneous so- 
lution. 

This paper presents a new approach to describe the 
widely scattered states for the synchronized mode of 
types (i) and (ii) where the multilane car correlations 
are fundamental in understanding the origin of conges- 
tion. The notion 'synchronized' traffic seems to embrace 
distinguished modes of congested multilane traffic flow. 
So we just deal with another type of the synchronized 
mode. We think that the role of the multilane corre- 
lations in the "oscillating congested traffic" is either of 
minor importance or is reduced to prevent the traffic flow 
going into the jammed "stop-and-go way" state. In the 



next section we will try to show that there is another 
type of synchronized mode further called 'light' synchro- 
nized traffic that meets the aforementioned Kerner's hy- 
potheses. Namely, in this mode the traffic flow dynamics 
contains fragments at which the vehicle velocity keeps 
its value approximately unchanged inside relatively long 
time intervals. The transitions between these fragments 
proceed through sharp jumps and the particular values 
of such a quasi-stationary vehicle velocity vary inside a 
wide interval. The subsequent sections will present the 
corresponding model developed from our approach for- 
mulated in papers 
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It should be noted that in the 
presented model these quasi-stationary states are not sta- 
tionary at all, at least in the standard physical meaning, 
so we prefer to call them long-lived states of synchronized 
traffic. 



II. SINGLE- VEHICLE DATA AND TWO TYPES 
OF THE OBSERVED SYNCHRONIZED TRAFFIC 



In this section we analyze single-vehicle data col- 
lected by loop detectors at the German highway Al near 
Cologne between 6-th and 17-th June 1996. Figure || 
illustrates the analyzed section of the highway and the 
position of the detectors. A detailed description of this 
data set and the detection devices has been presented 
in Jl6| |. Here, we only briefly recall their main features. 
The two sets of detectors Dl and D4 were placed nearby 
the busy intersection between the highway Al and the 
highway A 57 (Koln-Nord), while the detectors D2 and 
D3 are located close to the on-ramps and exit-ramps of 
the junction K61n-L6venich. In between there is another 
junction (AS K61n-Bocklemiind) but with a rather low 
usage. The most pronounced congested traffic was ob- 
served at the detector Dl where the number of lanes is re- 
duced from three to two for cars passing the intersection 
towards K61n-L6venich. Therefore, the data collected by 
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the selected detector were mainly studied in |16|]. The 
congested traffic observed near the detector Dl is char- 
acterized by a sharp fall in the mean velocity persisting 
within several hours. In connection with the following we 
prefer to call it heavy synchronized mode or "oscillating 
congested traffic" . 

The detectors D2, D3, and D4 also recorded congested 
traffic that, however, seems to be of another nature that 
will be called light synchronized mode in the following. 
Exactly in the vicinity of the detectors D2, D3, and D4 
traffic flow demonstrated the behavior meeting Kerner's 
hypotheses in the sense discussed at the end of the pre- 
vious section. The present paper does not pretend to 
a thorough analysis of these data, we only demonstrate 
the fact that the hypothetical behavior of traffic flow has 
been recorded directly and is likely not to be a misinter- 
pretation of non-stationary inhomogeneous vehicle struc- 
tures moving on the highway. 

For each detector the data set is mainly composed of 
the sequence of numbers {ti, m, U, fci} showing that i-th 
car of type fcj (passenger car, truck, truck trailer etc.) 
passed the given detector at time ti moving on lane U 
with speed Uj. It should be pointed out that the passing 
times {ti} were measured with second's accuracy. Below, 
the time variations of qi (*, 0) and vi (t, 0) will be ana- 
lyzed. Here, qi(t,Q) is the traffic flow rate and vi(t, 0) 
the vehicle velocity on lane I at time t (measured in sec- 
onds) averaged over the time interval 20 according to: 



ffl (t,0) = 5^u> e (t-ti), 

ti 

qi(t,Q) ^ 

The weight coefficients w%{£) are given by: 



(1) 
(2) 



we(t) 



Zexp (-|t|/9) for |*| < 20, 
for 1*1 > 20, 



(3) 



subject to the normalization condition 

oo 

E »eW = l, (4) 

t— — OO 

For > 1, the factor Z is given by Z ~ [20(1 - e" 2 )]" 1 . 
According to this expression the characteristic time scale 
of averaging with the weight u>e(*) is 20, that is why we 
will refer to the time interval of averaging the data set 
keeping in mind the value 20 rather than 0. 

As has been mentioned above, heavy congested traffic 
is a typical case in the vicinity of the detector Dl during 
the morning hours of workdays and is accompanied with 
a substantial fall in the mean vehicle velocity. The other 
detectors measured solely light congested traffic when the 
mean velocity, on the average, does not drop essentially 
but exhibits strong fluctuations only. To illustrate these 
types of traffic flow we have chosen, by way of example, 



the time series of velocity and flow collected by the de- 
tectors Dl and D2 on Th 6 June 1996 and Sa 15 June 
1996, respectively. Figure || visualizes these data. The 
row 1 left and row 2 left windows show time variations 
in the vehicle velocities at different lanes averaged over 
1 minute time interval according to Eq. (|J). The visu- 
alization time intervals (from 9:30 to 12:30 a.m. for D2 
and from 9:00 to 12:00 a.m. for Dl) span all the time 
during which the congested flow was recorded on these 
days. In spite of the clearly visible difference in the time 
patterns at the detectors D2 and Dl the vehicle veloc- 
ities at the left and right lanes are strongly correlated. 
So both the observed modes of congested traffic can be 
categorized as the synchronized traffic flow. It should be 
noted that the mean vehicle speed at the off-ramp lane 
near the detector Dl (Fig. ||) was at least twice as high 
as the mean speeds at the inner-lanes on 6 June 1996 and 
did not exhibit strong fluctuations during the given time 
interval. It seems that the traffic flow on this off-ramp 
was "free" in comparison with the main stream and did 
not affect it essentially. 

The row 1 right and row 2 right windows demonstrate 
time variations in the vehicle velocity and the traffic flow 
rate within two-lane 4-minute averaging. The depicted 
time interval has been extended to make possible com- 
parison of the congested state with the free flow. We see 
that the flow rate at the detector D2 did not practically 
show any visible dependence on the traffic state, whereas 
in the vicinity of the detector Dl it exhibits strong vari- 
ations rather than a considerable fall. The latter feature 
prompts us that the congestion observed near the detec- 
tor Dl is not the developed "stop-and-go wave" pattern. 
This statement, however, will be justified more carefully 
below. 

Let us now demonstrate that the two time series 
recorded at Dl and D2 present actually two synchronized 
modes different in nature. To substantiate this statement 
we depicted in the row 3 left and right windows the same 
time series of the vehicle velocity but on smaller scales 
(from 10:00 to 11:00 a.m.) and averaged only over 30- 
second interval. First, we see that the multilane corre- 
lations hold also on these scales. Second, the difference 
between the two time series becomes pronounced. The 
time series collected by D2 (the row 3 left window) con- 
tains several quasi-plateaus in the time variations of the 
vehicle velocities (singled out in this window and also in 
the row 4 left one by dotted ellipses). In contrast, the 
velocity time series obtained by Dl seems to be no more 
than a collection of spikes varying in amplitude. To ar- 
gue for the statement about the different nature of these 
congestions the row 4 left and right windows compare 
these 30-second averaging time series of the vehicle ve- 
locity and flow rate simultaneously for the left lanes. For 
the detector D2 we meet a visible distinction between the 
time pattern of vehicle velocity and that of the flow rate. 
For the detector Dl they are quite similar in appearance. 

To make this difference more pronounced we have an- 
alyzed the variance for the time series of the vehicle ve- 
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FIG. 3: Time variations in the vehicle velocity and flow rate observed at the detectors Dl and D2 on Th 6 June 1996 and Sa 15 
June 1996, respectively. This figure visualizes the characteristic properties of two limit types of the synchronized mode, lightly 
congested traffic (D2) and heavily (Dl) one. Red color is chosen to lighten the comparison of the vehicle velocity dynamics 
at the left and right lanes, the blue curves are the flow rate dynamics. Dotted ellipses single out possible fragments of time 
variations in the vehicle velocity that can be referred to as the long-lived states of the synchronized mode. In the text particular 
windows are labeled with the row number counted from the top and the side (left or right). 
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lane per two seconds. Keeping in mind that the detectors 
measured the times of car passing within seconds we may 
treat each record of a car having passed the detector as 
a random event with the probability p ~ 0.5. Let us 
consider these events mutually independent and describe 
them with a sequence {7(2)} of random numbers such 
that take values equal to either 1 or with the probability 
p and 1 — p, respective, i.e. 

(7(0) =P and (7(ih(i'))=pSu'+p 2 (l-6u'), (6) 

where Su> is Kronecker's symbol. In these terms Eq. (|l|) 
can be rewritten as 



FIG. 4: The relative deviations of the vehicle velocity and flow 
rate vs the averaging interval for the time series collected by 
the detectors Dl and D2 at the left lane on 6 June 1996 and 
15 June 1996, respectively, between 10 and 11 a.m. 



locity and flow rate collected by the detectors Dl and D2 
at the left lane during the observation time from 10:00 
to 11:00 a.m. Fig. [| presents the corresponding relative 
deviations 
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depending on the averaging interval 20. For the Dl time 
series we meet at first a weak dependence of the relative 
deviations on the averaging scale 20 for both the vehicle 
velocity and flow rate which, then, decrease fast keeping 
approximately their ratio fixed, i.e. 



Sv(Q) 
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However, the question as to whether the latter estimate 
has a certain physical meaning or just is a simple coinci- 
dence requires an individual analysis. In contrast, the D2 
time series demonstrates another behavior. The relative 
deviation of the vehicle velocity exhibits a weak decrease 
approximately from 22% to 13% as the averaging interval 
grows 16-fold from 30 seconds to 8 minutes. Conversely, 
the relative deviation of the flow rate is scaled practically 
as (straight solid line in Fig. ^) 



Sq(Q) 



(q) 



12% x 



D2 



0.5 min 
26 



(5) 



Expression (|5|) enables us to suppose that the fluctuations 
of the flow rate observed at the detector D2 during one 
hour from 10 to 11 a.m. are due to a certain white noise. 
Moreover, the mean flow rate (q)\- D2 left lane during the 
given time interval was about 1800 veh/h/1, so on average 
the detector D2 recorded one car passing it on the left 



Whence we immediately get the following expression 
for the mean relative deviation of the random variable 



(([ gi ft,e)-( g ^,9))] 2 )) 

(9l(i.6)) 



2p 



5T (7) 



In obtaining Eq. (0) we have taken into account the nor- 
malization condition (^) and for the sake of simplicity 
extended the exponential dependence (||) of the weight 
coefficients we,(t) to all the values of its argument t. 
It should be recalled that the value appearing in (||) 
matches the time scale measured in seconds. For p = 0.5 
and the time scale measured in minutes Eq. ([?]) directly 
gives us expression with the replacement of 12% by 
13%. 

Thereby, the time variations in the flow rate recorded 
by the detector D2 are likely to be due to the discrete- 
ness of traffic flow bearing no correlations at all. In no 
case such fluctuations should be ascribed to the state of 
traffic flow. Thus dealing with the congested traffic in 
the vicinity of the detector D2 we actually meet the traf- 
fic flow state where the vehicle velocity exhibits strong 
long-time fluctuations with certain time plateaus and the 
traffic rate, in contrast, has to be regarded as a constant 
value. These plateaus are separated by sharp jumps as 
illustrated in Fig. ||, the row 3 left and row 4 left win- 
dows with respect the plateaus singled out by the first 
two dotted ellipses. In the row 3 left windows it is also 
clearly visible that the plateau position with respect to 
the u-axis varies substantially. It should be noted that 
Neubert et al. jl6) also have found strong long-time cor- 
relations in the vehicle velocity and their absence in the 
headway distribution but with respect to the Dl time se- 
ries. They also have obtained the flow rate independence 
from the vehicle density for the synchronized mode but 
after averaging over all the Dl time series. Besides, there 
are cellular automata models predicting also a constant 
value of the flow rate inside a certain interval of the ve- 
hicle density which is due to local highway defects (see 
Ref. [1| [47|, g |49|, |0| and the review by Chowdhury et 
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FIG. 5: The fundamental diagrams showing the possible states of traffic flow in the pq-plane where for a chosen averaging 
interval 20 the vehicle density p(t, Q) was calculated as the ratio of the current flow rate q(t, O) to vehicle velocity v(t, 0), i.e. 
p = q/v. The analyzed time series (visualized in Fig. |^) and the corresponding colors are depicted at the window tops. Red 
color is used for the traffic flow states recorded within the time interval [10:00,11:00] a.m. (the corresponding time series are 
shown in Fig. || the row 3, 4, left, right windows). Rows 1 and 2 illustrate evolution of the fundamental diagrams depending 
on the averaging interval, whereas the windows in the third row show the [10:00,11:00] series in detail individually. In these 
windows the red dots are time series points connected by solid black lines according their order in the time series. Blue lines 
are drawn to guide the eyes to follow the quasi-free traffic states specified by the relationship q — vp between the flow rate q 
and the vehicle density p where the coefficient v is a certain constant. 
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al. ||). However, in the case under consideration a con- 
stant value should be ascribed to the current flow rate. 
The latter enables us also to suppose that a possible spa- 
tial pattern of the vehicle distribution was fixed on the 
highway during the observation time. So the D2 time 
series of the vehicle velocity is likely to visualize the real 
dynamics of traffic flow caused by its internal local prop- 
erties rather than by the pinch effects. In this connection 
it should noted that Lee et al. j3^, [h| and Helbing et al. 
jllj preferred to place the congested traffic states match- 
ing spatial vehicle patterns fixed as a whole on highways 
into an individual group. 

In order to demonstrate that the analyzed traffic states 
belong in fact to the synchronized mode we present the 
fundamental diagram (Fig. ||) reflecting the traffic flow 
states on the pg-plane. For a chosen averaging inter- 
val 29 the current values of the flow rate q(t, 0) and 
the vehicle velocity v(t, 6) were calculated according to 
Eqs. (|IJ), (|^), and then the effective vehicle density p 
was calculated as p(t,&) = q(t,Q)/v(t,Q). The row 1 
and row 2 windows exhibit the fundamental diagrams 
for the traffic flow on the left lane near the detectors D2 
and Dl, respectively, depending on the averaging interval 
26. Red color has been used to make visible the states 
of traffic flow recorded within the time from 10 a.m. to 
11 a.m. and thus to separate the free flow states from the 
those of the congested traffic. We see that, in fact, the 
congested traffic flow under consideration possesses the 
widely scattered state domain and so can be referred to 
as the synchronized mode. It is most clear for the time 
series averaged over 1-minute interval. On one hand, the 
widely scattered states cover a sufficiently large domain 
on the pg-plane. On the other hand, for the D2 time se- 
ries the gap between this domain and the free flow region 
becomes quite visible. For the Dl time series a certain 
"tail" of the widely scattered state domain in the vicinity 
of the origin (p = 0, q = 0) (see the row 2 left window) 
typical for the "stop-and-go way" pattern [jl6| disappears 
on these scales. So the "stop-and-go way" traffic flow was 
not, at least, the dominant phase state. 

The windows in row 3 show the fundamental diagram 
for the traffic flow states recorded during the time from 10 
to 11 a.m. These windows actually exhibit the structure 
of the widely scattered state domains visualized in the 
given time series. Roughly speaking, both of them com- 
prise a large amount of straight lines describing quasi- free 
flow states characterized by the relationship q = pv with 
a constant coefficient v and the continuous transitions be- 
tween them. The particular magnitudes of v, however, do 
not coincide with the free flow velocity being about 100 
km/h (see Fig. [|, rows 1,2 right windows) but vary inside 
a wide interval from 35 km/h to 90 km/h for the D2 series 
and from 5 km/h to 45 km/h for the Dl series. Concern- 
ing the lower boundary of the latter interval, however, it 
should be pointed out that the result can depend on the 
way of calculating the vehicle density |L6). For the D2 
time series this feature remains also after averaging over 
1-minute interval, whereas for the Dl series it disappears 



and only the transition lines are visible. 

The row 1 and row 2 windows again raise the ques- 
tion about the fundamental diagram which has been dis- 
cussed in many papers (see, e.g., reviews in partic- 
ular, which averaging interval should be used. As can be 
seen in these windows the scattered state domain shrank 
considerably on scales about 8 minutes. Besides, keep- 
ing in mind the drawn conclusion about the nature of 
the flow rate fluctuations recorded by the detector D2 
we can suppose that the widely scattered states of the 
light synchronized mode are due to the traffic flow dis- 
creteness. The vehicle velocity, however, exhibit strong 
fluctuations comprising many long-lived states and the 
continuous jump transitions between them. 



III. SYNCHRONIZED MODE AS A 
GLASS-LIKE PHASE STATE 

When vehicles move on a multilanc highway without 
changing the lanes they interact practically with the 
nearest neighbors ahead only. Therefore, there should 
be no internal correlations in the vehicle flow at differ- 
ent lanes. In particular, the drivers that would prefer to 
move faster than the statistically mean driver will bunch 
up forming the platoons headed by a relatively slower 
vehicle. When the cars begin to change lanes for over- 
taking slow vehicles the car ensembles at different lanes 
will affect one another. The cause of this interaction is 
due to the fact that a car occupies during a lane change 
manoeuvre two lanes simultaneously, affecting the cars 
moving behind it at both lanes. Figure ||(6) illustrates 
this interaction for cars 1 and 2 through car 4 changing 
the lanes. The drivers of both cars 1 and 2 have to regard 
car 4 as the nearest neighbor and, so, their motion will 
be correlated during the given manoeuvres. In the same 
way car 1 is affected by car 3 because the motion of car 
4 directly depends on the behavior of car 3. The more 
frequently lane changing is performed, the more corre- 
lated traffic flow is on a multilane highway. Therefore, 
to characterize traffic flow on multilanc highways it is 
reasonable to introduce an additional state variable, the 
order parameter h \ 



45 1 . This variable is the mean 



density of such car triplets normalized to the maximum 
possible for the given highway and, so, can play the role 
of a measure of the multilane correlations in the vehicle 
flow. 

On the other hand, the order parameter h introduced 
in this way can be regarded as a measure of the vehi- 
cle arrangement regularity. Let us discuss this question 
in detail for the free flow and synchronized mode indi- 
vidually. In the free flow the feasibility of overtaking 
makes the vehicle arrangement more regular because of 
platoon dissipation. So, as the order parameter h grows 
the free traffic becomes more regular. Nevertheless, in 
this case the density of the car multilane triplets remains 
relatively low, /i<l, and the vehicle ensembles should 
exhibit weak correlations. Figure [7] illustrates the F — > S 
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order parameter h : for (a) and (c) h <C 1 ; for (b) h ~ 1 
(a) low car density, p«p (c) high car density, p ~ p 

^ od-» id— od 

CD) HO* HO* E2 



(6) car density near the transition "free flow -** synchronized mode", p ~ p c 




FIG. 6: Schematic illustration of the car arrangement in the 
various phases of traffic flow and the multilane vehicle inter- 
action caused by cars changing the lanes. 



(a) Free flow: car density p « p c , order parameter kf< h s 




(b) Synchronized mode: car density p » p , order parameter h s > hf 




FIG. 7: Schematic illustration of the alteration in the vehicle 
arrangement near the "free flow <-» synchronized mode" phase 
transition. 



transition. As the car density grows in free flow, the 
"fast" drivers that at first overtake slow vehicles indi- 
vidually begin to gather into platoons headed by more 
"slow" cars among them but, nevertheless, moving faster 
than the statistically mean vehicle, see Fig. |(a). The 
platoons are formed by drivers preferring to move as fast 
as possible keeping short headway distance without lane 
changing. Such a state of the traffic flow should be suf- 
ficiently inhomogeneous and is the reason for the occur- 
rence of high-flow states in free traffic. Therefore, even 
at a sufficiently high car density the free flow should be 
characterized by weak multilane correlations and not too 
big values of the order parameter hf. The structure of 
these platoons is also inhomogeneous: they comprise cars 
whose drivers would prefer to move at different headways 
(for a fixed velocity) under comfortable conditions. So, 
when the density of vehicles attains sufficiently high val- 
ues and their mean velocity decreases remarkably with 
respect to the velocity on the empty highway some of 
the "fast" drivers can decide that there is no reason to 
move so slowly at such short headway requiring strain. 



Then they can either overtake the car heading the cur- 
rent platoon by changing lanes individually or leave the 
platoon and take vacant places, compare Fig. 0(a). The 
former has to increase the multilane correlations and, in 
part, to decrease the mean vehicle velocity because the 
other drivers should give place for this manoeuvres in suf- 
ficiently dense traffic flow. The latter also will decrease 
the mean vehicle velocity because these places are vacant 
from the standpoint of sufficiently "fast" drivers only but 
not from the point of view of the statistically mean ones 
preferring to keep longer headway in comparison with 
the platoon headway. That means, that the statistically 
mean drivers have to decelerate, decreasing the mean ve- 
hicle velocity. The two manoeuvre types make the traf- 
fic flow more homogeneous by dissipating the platoons 
and smoothing the headway distribution, see Fig. 0(6). 
Besides, the single-vehicle data [[l6) show that the syn- 
chronized mode is characterized by long-distant correla- 
tions in the vehicle velocities. The headway fluctuations 
are correlated only on small scales. These findings jus- 
tify the assumptions of the synchronized mode being a 
more homogeneous state than the free flow. We think 
J45[ that the given scenario describes the synchronized 
mode formation which must be characterized by a large 
value of the order parameter, h s > hf, and a lower veloc- 
ity in comparison with the free flow at the same vehicle 
density. It should be noted that this mechanism of the 
synchronized mode emergence is related to the hypoth- 
esis by Kerner ||, ||l| that the synchronized mode is 
caused by a "Z"-like form of the overtaking probability 
depending on the vehicle density. 

Keeping in mind this scenario we have assumed the 
mean velocity v of multilane traffic flow to be determined 
by both the vehicle density p and the order parameter 
h, namely, v = p), where the function §{h,p) is 
regarded as a phenomenological relationship known be- 
forehand. Its general properties have been discussed in 
Ref. HJ, |^|, in particular, it has been assumed to be a de- 
creasing function with respect to h for the synchronized 
mode. Previously we have written the governing equation 
for the order parameter completely describing the traffic 
state evolution in terms of time variations in the mean 
velocity v and the order parameter h, the vehicle density 
is treated as a fixed value. In this way we have obtained a 
simple and natural explanation of the observed sequence 
of phase transitions F — > S — > J, with each of them be- 
ing of the first order. In particular, this model predicts 
that the order parameter h exhibits a sharp jump from 
the value hp(p c ) -C 1 corresponding to the free flow up 
to the order parameter hsific) ~ 1 °f the synchronized 
mode when the vehicle density p exceeds a certain critical 
value p c . The further evolution of the synchronized mode 
as the vehicle density grows is given via the dependence 
v = $[h S (p),p\. 

However, to describe the widely scattered state region 
of the synchronized mode we have to regard the multi- 
lane car interaction into more detail. The matter is that 
for a car to be able to leave a platoon or to change a 
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Synchronized mode: car arrangement enabling (a) or hindering (b) overtaking 



^> BCD' — ► 




FIG. 8: Schematic illustration of the car arrangement at the 
neighboring lane in the synchronized mode that enables over- 
taking (a) or hinders it (6). 



lane, the local vehicle arrangement at the neighboring 
lanes should be of a special form (Fig. ||). So the car 
rearrangement essentially depends also on the particu- 
lar details of the neighboring car configuration exhibiting 
substantial fluctuations. These fluctuations, however, are 
of another nature than those in the free flow. The lat- 
ter can be treated as a white noise and slightly blurs the 
dependence v = d(h,p) for the free flow (Fig. [l]). In the 
synchronized mode, by contrast, attaining the optimal 
conditions of driving, including also overtaking slow ve- 
hicles, can be frustrated for a certain time as illustrated 
in Fig. H For car 1 to be able to overtake car 2 the neigh- 
boring car 3 should provide a room for this manoeuvre. 
Otherwise the driver of car 1 has to wait and the lo- 
cal car arrangement will not vary substantially. In other 
words, changes in the particular realizations of the local 
car arrangement can be frozen for a certain time although 
the globally optimal car configuration is not attained at 
the current time moment. Due to this self-freezing ef- 
fect the synchronized mode can comprise a great amount 
of locally metastable states and correspond to a certain 
two-dimensional region in the pg-plane rather than to a 
line q = fl(p)p. This feature seems to be similar to that 
met in physical media with local order, for example, in 
glasses where phase transitions are characterized by wide 
range of controlling parameters (temperature, pressure, 
etc.) rather than their fixed values (see, e.g., Q). 

The synchronized mode consists of clusters of some ten 
cars moving along the road as a whole jl6| which will be 
called for convenience the fundamental clusters of the 
synchronized mode. Thereby, the introduced order pa- 
rameter h is actually an averaged characteristics of such 
a cluster as whole and, so, it does not allow for local 
fluctuations in the car arrangement inside this cluster. 
To describe the latter another state variable, the order 
parameter a has to be introduced. Its physical meaning 
is to allow for the fact that the optimal and comfort- 
able way of driving which the drivers try to attain indi- 
vidually depends essentially on the current realization of 
the car arrangement inside the given fundamental clus- 
ter. In other words, we assume that the order parameter 
/io(a,p) matching this optimal way of driving is a func- 



tion of the vehicle density p as well as the order parameter 
a. Without lost of generality in such a phcnomenologi- 
cal approach we may relate ho(a = 0,p) to the global 
optimum of the driving manner attained after averaging 
over all the realizations of the car configuration, i.e. set 
ho(a — 0, p) — hs{p). Besides, it is reasonable to relate 

variations of the order parameter a of amplitude about 

l li 

unity, (a 2 ) ~ 1, with all the possible realizations of 
the car configuration. 

To complete the description of the long-lived state con- 
tinuum of the synchronized mode we should specify the 
evolution of the order parameter a. Since the order pa- 
rameter a allows for microscopic details of the funda- 
mental cluster structure its fluctuations will be treated 
as a random noise whose amplitude depend on the ve- 
hicle density only. In contrast, the rate da/dt of time 
variations in the order parameter a has to be affected 
substantially by the current value of the order parameter 
h. In fact, as the order parameter h tends to the local op- 
timum value /io(a, p) for the given a the rate da/dt should 
be depressed because all the drivers forming the funda- 
mental cluster prefer to wait until a more comfortable car 
configuration arises therefore inhibiting the evolution of 
the fundamental cluster structure. 

Now, the model can be formulated. However, before 
doing this it is worth noting that the present paper does 
not describes the formation and evolution of the synchro- 
nized mode inside the free flow. We construct a contin- 
uum model for the long-lived states. So we may confine 
our consideration solely to the local properties ignoring 
the interaction between fragments of the synchronized 
traffic flow separated in space. The latter is worthy of 
individual and detailed investigations. 

IV. THE DYNAMICAL TRAP MODEL 

Keeping in mind the aforesaid and following Q, p5[ , 
we specify the local state of traffic flow on a multilane 
highway by the system of four phase variables, the mean 
velocity v, the density p of vehicles, and two order pa- 
rameters h and a. The parameter h describes the effect of 
multilane car interaction on the individual driver behav- 
ior, so, the mean velocity v = $(h, p) is a direct function 
of the order parameter h and the vehicle density p given 
beforehand and meeting general properties described in 
^5) . The other order parameter a allows for fluctua- 
tions in the car arrangement inside a fundamental cluster 
of the synchronized mode and affects the relaxation pro- 
cess of the order parameter h, namely, we write 

T^ = -<f>(h,a,p) +eV^(*)- (8) 

Here, r is the average time drivers require to come to 
the decision to begin or to stop overtaking manoeuvres, 
the function <p(h, a, p) describes the regular behavior of 
individual drivers depending on the traffic state as well 
as on the particular details of the local car arrangement. 



11 



Finally, the term ey/r^(t), a stochastic Langevin force, 
allows for a random component in the driver evaluation 
of the current situation on the road: 



(9) 



The factor e quantifies the size of the noise. As dis- 
cussed in the previous section it is reasonable to nor- 
malize the order parameters h and a to their maximum 
values, which enables us to accept that in the free flow 
the averaged value of the order parameter h p [p] is much 
less than unity, hp[p] -C 1, whereas in the synchronized 
mode it is about one, hs[p] ^ 1, and to assume that the 
order parameter a varies around zero, i.e. (a 2 ) ~ 1. 

In order to proceed, recall briefly the key features of 
the synchronized mode. When the vehicle density ex- 
ceeds a certain critical value p c the free flow becomes 
unstable and the synchronized modes arises on the road 
through a first order transition. Our previous model de- 
scribes this F — > S transition as a sharp jump in the 
order parameter h from hp[p] to hs[p] at a fixed value of 
the vehicle density p which should be accompanied with 
stepwise jumps in the mean velocity 8v — vp[p] — vs[p] 
and the traffic flow rate Sq — p ■ Sv. These jumps stem 
from the dependence of the steady state order param- 
eter h st [p] on the vehicle density p taking the multi- 
valued form for p > p c . In the given case the points 
(p,qp[p\) and (p, qs[p\) on the pq-plane have to corre- 
spond to well determined states of traffic flow. However, 
the available empirical data (illustrated in Fig. justify 
the latter assumption with respect to the free flow only, 
whereas the synchronized mode stretches over a whole re- 
gion with, may be, blurred boundaries qjj [p] and qg [p]- 
Besides, this region turns out to be sufficiently thick, 

\p\ ~ Qs '\fi\) ~ 1f[Pc] and adjacent to the free flow 



branch, ( qp[p] — c?g [p]) <C qp[p c ]- We intend to ex- 
plain the existence of the synchronized mode continuum 
by the effect of random fluctuations in the vehicle veloc- 
ity v. However, it is not sufficient to introduce additional 
stochastic Langevin forces into governing equations, as 
done with respect to equation (B). Indeed, in this case 
such a stochastic force quantifies the random component 
of the individual driver behavior, so, it should demon- 
strate similar features in both the free flow and synchro- 
nized mode. Therefore, in particular, we consider the 
amplitude e of the introduced stochastic force as well as 
its correlation function (||) to be independent of traffic 
flow state. The free flow matches a line on the pq-pl&ne, 
the free flow branch, which is slightly blurred due to ran- 
dom fluctuations. Besides, it exhibits no long-time cor- 
relation. Thereby, an additional Langevin force cannot 
explain the wide quasi-equilibrium state continuum of 
the synchronized mode. Moreover, these Langevin forces 
cannot be strong, so, we may set e -C 1 in equation 

Furthermore, we assume that in the synchronized 
mode where the multilane car interaction is essential the 
order parameter h is substantially affected by the order 
parameter a quantifying the random fluctuations in the 



local car arrangement. In this case the possibility of a 
self-freezing effect discussed above is responsible for the 
emergence of the long-lived state continuum. Therefore, 
the next step should specify the function 4>(h, a, p) in such 
a way that the resulting steady-state dependence h st (p, a) 
on the order parameter a is strong. Since we deal solely 
with the synchronized mode itself we can write down 



cf>(h,a,p) =h - h (a,p) 



(10) 



without lost of generality. Here the function ho(a, p) can 
be approximated by the expression: 



h (a,p) = h s (p) - Aa 2 



(11) 



which assumes that a zero value of the order parame- 
ter a corresponds to the most comfortable conditions for 
driving, to the highest order of the traffic flow, and to 
the slowest vehicle motion. We note that in the previous 
model 1 44, E5j the value hs{p) characterizes the station- 



ary state of the synchronized mode at the given vehicle 
density p. The deviation of the car arrangement from the 
optimal configuration destroys the order of traffic flow, 
bringing the synchronized mode closer to the free flow 
in properties. Therefore, we may set A approximately 
equal to unity, A ~ 1, and to regard the traffic state 
with h w ft,o(l,p) as practically the free flow or, more 
rigorously, the state directly adjacent to the free flow. 

Now let us specify the governing equation for the or- 
der parameter a. Local variations of the car arrangement 
seem to be determined solely by the vehicle density. In- 
deed, the order parameter h is actually produced by av- 
eraging, first, over many particular realizations of the 
fundamental cluster structure similar in properties with 
respect to the possibility of the lane changing and, sec- 
ond, over all the cars making up the fundamental cluster. 
So it is sufficiently rough characteristics of the traffic flow 
state and the dynamics of the order parameter a can be 
treated as an independent of h random process of unit 
amplitude and induced by the white noise except for the 
rate Q(h, a) at which the parameter a responds to the 
white noise. Namely, we write 

T ^ = -n(h,a)a + n 1 / 2 (h,a)V?<;(t), (12) 
at 

where again the noise <;(t) meets the conditions 

(?(t)) = o, man) = o, (<r(*Mi')> - m n w 

and the factor D,(h, a) describes the self-freezing effect of 
fluctuations in the car arrangement. In particular, when 
the order parameter comes close to the local equilibrium, 
h — » ho(a,p), i.e. (j>(h,a,p) — > the rate Q,(h, a) should 
tend to zero. Therefore, by keeping only the leading term 
we get 



Q{h, a) 



4> 2 {h,a,p)/(j>l if \cj)(h,a,p)\ < 4> 
VL if \<p(h,a,p)\ > 0o 



(14) 
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Here, the constant coefficients c/)q ~ ilo ~ 1 because un- 
der the generaf conditions time variations of the behavior 
of individual drivers as well as of the car configuration 
of fundamental cluster should proceed essentially at the 
same rate. The form of the random force in equation Jl2] ) 
has been chosen so that the amplitude of fluctuations in 
variable a remains unchanged provided the factor Q(h, a) 
is constant. In other words, the proximity of the funda- 
mental cluster structure to a local equilibrium slows down 
the time variations in the order parameter a. In some 
sense the curve h = ho(a, p) on the a/i-plane (for a fixed 
value of p) forms certain traps for the path {a(t),h(t)}, 
which has given the model its name. If the order param- 
eter h had been constant, governing equation fll2| ) would 
describe the collapse of the random motion {a(t)}, with 
the Ito process corresponding to the maximum of the 
collapse intensity (see, e.g., Ref. |53|). Therefore we also 



ascribe to the stochastic equation (12) the Ito type, which 
completes the model formulation. Naturally, in the given 
model the vehicle density p is a fixed constant. Time vari- 
ations in the order parameter h(t) can be mapped onto 
the evolution of the traffic flow rate q(t) = pi3(h(t), p). 
Therefore, in the following it is sufficient to regard the 
order parameters a and h only. 

Here we present only the preliminary investigation of 
the stated model, namely, the results obtained by its nu- 
merical simulation, which demonstrate that this model 
in fact leads the desired long-lived state continuum. For 
the sake of continuance in numerical simulations we have 
converted to the dimensionless time t — > t/r, introduced 
the variable 77 = (hs{p) — h)/A, renormalized the con- 
stants e — > e/A and 4>o — ¥ 0oA, set = 1. This leads 
to: 



da 



em 



(15) 



— = -lW(r/,a)a + fto /2 w 1/2 (77, ak(i) (16) 

where the noise and meets the conditions @ 
and (|l3|). Finally, the function 07(77, a) may be specified 
as: 



w(r], a) 



(ry-a 2 ) 2 /^ if 



1 



if 



Again, the stochastic equation ( [Tq ) is assumed to be of 
the Ito type. 

Let us now discuss the results obtained by simulating 
numerically the system of Eqs. ([t5]), (O) for Oo = 1 , 
e = 0.1, and cf>o = 0.5. First of all, Fig.^jj illustrates the 
evolution of the order parameters 77 and a when there is 
no self- freezing. In this case the order parameter a ex- 
hibits the standard random pattern of Brownian move- 
ment inside a region of unit width. We see a collection of 
practically independent spikes of unit width. A similar 
pattern (see Fig. [|(a)) is demonstrated by the dynamics 
of the order parameter 77 provided the interaction of the 
parameters 77 and a have been ignored, i.e. by replacing 
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by 77. Naturally, in this case the synchronized 



FIG. 9: The time pattern of the order parameters r\ and a if 
the self-freezing effect is suppressed. The window (a) exhibits 
time variations in the order parameter 77 in the case where 
the influence of the order parameter a has been ignored at 
all. The window (b) presents the dynamics of 77 affected by a. 



mode matches a line on the pg-plane for e <C 1. The 
dynamics of the order parameter 77 affected by the vari- 
able a but without the reciprocal influence is shown in 
Fig. ^(6) . Again a collection of spikes can be seen, whose 
amplitude as well as width has increased tenfold. In other 
words, if we take into account the effect of the fluctua- 
tions in the car arrangement (the order parameter a) on 
the individual driver behavior (the order parameter h) 
then it is possible to explain an essential blurring of the 
synchronized mode state. 

The dynamics changes dramatically for the full prob- 
lem, see Fig. [l(]. The time pattern takes a form corre- 
sponding to the long-lived state continuum. When the 
point {a(t), T](t)} representing the current state of the 
synchronized mode wanders on the a?7-plane and reaches 
the curve 77 = o 2 at any point it will be trapped for a cer- 
tain time until it finally escapes from the trap due to the 
noise e£(t). After that the system again wanders in the 
a?7-planc during a time interval about unity before being 
trapped for the next time. Since the characteristic dura- 
tion of the trapping is much longer than unity the pattern 
looks like a certain collection of local metastable states of 
the synchronized mode. However, such prolonged stays 
of the system are not metastable in the rigorous mean- 
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FIG. 10: The time pattern of the order parameters rj and a 
when the self-freezing effect is substantial. 



ing and we preferred to call them simply the long-lived 
states. Since each point of the curve r\ = a 2 is a trap, the 
long-lived states make up a certain continuum. 



V. REMARKS ABOUT GOVERNING 
EQUATIONS FOR LIGHT SYNCHRONIZED 
TRAFFIC 

In writing equations ( p"5j ) and ( |l6| ) we actually have 
considered the state of the traffic flow homogeneous. 
Thereby, the interactions between different fragments 
separated in space have been neglected. In order to an- 
alyze a more realistic case of the evolution of traffic flow 
with the long-lived states we have to deal with the fields 
p(x, t), rj(x, i), and a(x, t) describing the vehicle distribu- 
tion along the highway, the order of the car arrangement 
as well as local fluctuations in the car arrangement at dif- 
ferent highway parts. It should be noted that the given 
model has been proposed for the light synchronized mode 
where the driver's delay to the headway variations seems 
to be of minor importance. This enables us to assume 
the local value of the vehicle velocity practically equal 
to the optimal velocity ${h, p) for current values of the 
vehicle density p and the order parameter h. 

In this section we only touch a possible way of con- 
structing the distributed model for traffic flow with long- 
lived states. For inhomogeneous traffic flow Eqs. (15), 
( |l6| ) actually describe the process in the frame locally 



"attached" to the traffic flow, i.e. moving at the same 
speed along the road. So, in order to return to the phys- 
ical coordinate system we, first, should replace the time 
derivative by the material derivative, i.e. 

— => d t + u(r), p)d x , 

where we again use the dimensionless form of time t in 
the term u(j], p) = T-d[{hs{p) — A?y), p\. 

Second, as was discussed in Sec. [nj the order param- 
eter h characterizes the arrangement of the fundamental 
car cluster of synchronized traffic flow. So it cannot ex- 
hibit considerable variations on spatial scales about the 
characteristics size I of this cluster. So the right-hand 
side of equation (|l5|) has to be modified to allow for this 
effect. Following our previous paper |45| we introduce in 
equation ( fl5| ) the term 

£ 2 dlh + 2- 1 / 2 £d x h, 

smoothing the field h on the scales I and taking into 
account the asymmetry in the vehicles interaction. Con- 
versely, the order parameter a allows for local variations 
of the car arrangement inside the fundamental car clus- 
ter, so we may regard the field a(x, t) as totally indepen- 
dent at different points along highway, at least, within 
such a mesoscopic description. 

In this way the system of equations Jl5| ) and ( |l6| ) can 
be extended as follows: 

d t ri + u(ri,p)d x ri = £ 2 d 2 h + 2~ 1 / 2 £d x h 

-(r,-a 2 )+eat,x) (17) 

difl, + u(r), p)d x a = —VL$m{ri,a)a 

+ n 1 /2 m 1 / 2 ( V ,a),(t,x),(18) 

where the terms £(t, x) and <r(t, a;) are random sources un- 
corrected in time and space. Naturally, these equations 
must be completed by the continuity equation: 



dtp + d x [u(r), p)p] = {Dr)d 2 x p , 



(19) 



where D is the effective diffusion coefficient. 

The proposed distributed model for the light synchro- 
nized traffic flow is worthy of an individual investigation. 
Here we have only demonstrated the feasibility of its con- 
struction. 



VI. CONCLUSION 

We proposed a novel notion of the widely scattered 
states in the synchronized traffic stimulated by Kerner's 
hypotheses ||, H 111' Hi about the existence of a 
multitude of metastable states in the homogeneous or 
"homogeneous-in-speed" synchronized traffic. 

First, we have analyzed the single- vehicle data col- 
lected at the German highway Al near Colonge from 6-th 
to 17-th of June 1996 in order to verify the presence of 
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the hypothetical multitude of metastable states. As a 
preliminary result we have found two types of synchro- 
nized mode observed at different detectors. One of them 
can be classified as the "oscillating congested traffic" |J 
and has been analyzed in detail in Eg ]. The behavior 
of the other synchronized mode differs essentially there- 
fore we preferred to classify it as an individual type called 
light synchronized traffic. Exactly this type of traffic flow 
possesses the hypothetical multitude of quasi-stationary 
states. Namely, it has been found that the vehicle veloc- 
ity dynamics consists of a collection of certain fragments 
inside which the vehicle velocity is approximately con- 
stant and the continuous transitions between them occur 
via sharp jumps. The particular values of the vehicle 
velocity in these fragments vary in a wide interval. In 
contrast, the traffic flow rate has to be regarded as a 
constant because, as we have shown, its fluctuations are 
mainly due to the traffic flow discreteness. In spite of 
the latter properties the fundamental diagram drawn for 
this traffic mode in the standard way contains the widely 
scattered state domain. Keeping in mind the obtained 
results we preferred to call this domain the long-lived 
state region, because it is made up of states that are not 
stationary at all, at least in the standard physical mean- 
ing. 

Second, we have proposed a simple mathematical 
model explaining this complex behavior of light synchro- 
nized traffic in terms of dynamical traps. Advancing the 
idea proposed in 



45 



we specify the state of traffic 
flow by four phase variables, the mean velocity v, the ve- 
hicle density p, and two order parameters h and a. The 
parameter h describes the effect of multilane correlations 
in the vehicle motion on the behavior of individual drivers 
originating from lane change manoeuvres. On the other 
hand, it characterizes the order of traffic flow, so the 
transition from the free flow to the synchronized mode is 
represented as a sharp jump in the value of the order pa- 
rameter h. We assume the existence of the relationship 
between the mean velocity v, the vehicle density p and 
the order parameter h, i.e. v = *&(h,p), where &(h, p) is 
a certain phenomenological function known beforehand. 



In these terms we can describe the local dynamics of the 
synchronized mode in terms of the evolution of the or- 
der parameter h, here the vehicle density p is treated 
as a fixed constant. The latter is justified because the 
synchronized mode is characterized by long-distance cor- 
relations, so its basic features should stem from the main 
properties of a certain fundamental cluster, the minimal 
fragment of car ensemble that moves as a whole and can 
be regarded as an atomic element of the synchronized 
mode. 

The dynamics of the order parameter h is governed 
by the vehicle density and in addition, by local fluctua- 
tions in the car arrangement affecting the feasibility of 
changing lanes for overtaking. The latter aspect is the 
key point of the mechanism responsible for the emer- 
gence of widely scattered synchronized states. Due to 
different realizations of the car arrangement inside the 
fundamental cluster there exist a continuum of locally 
quasi-equilibrium states of the order parameter h. So, 
when such a local quasi-equilibrium state is attained, the 
individual drivers prefer to wait for a more comfortable 
car configuration to continue overtaking. In this case 
also the rate of changes in the car arrangement slows 
down and, as a result, the current state of traffic flow is 
'frozen'. Exactly this effect of self-freezing gives rise to a 
continuum of long-lived states in the light synchronized 
traffic. 

We have written a particular mathematical description 
of self-freezing and justified these qualitative conclusions 
by solving the stated model numerically. 
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